Spectroscopy of Superradiant scattering from an array of Bose-Einstein Condensates 
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We theoretically study the superradiant gain and the direction of this coherent radiant for an 
array of Bose-Einstein condensates in an optical lattice. We find that the density grating is formed 
to amplify the scattering light within the phase match condition. The scattering spectroscopy in 
the momentum space can provide a method for measuring the overlap of wavefunction between the 
neighboring sites, which is related to their inner-site and inter-site coherence. 



PACS numbers: 



I. INTRODUCTION 

Superradiance from a Bose-Einstein Condensation 
(BEC) offers the possibihty to study the novel physics 
associated with cooperative scattering of light in ultra- 
cold atomic systems. A series of experiments n], l3. |3, 
i, H i, 0, 1 and related theories [i,ll3, [IH, EllI llil 
have sparked the related interest in the quantum infor- 
mation , collective instability [l^, [l^ , high precision 
measurement 17 1 , and coherent atom optics Il8| . 

In a typical BEC superradiant experiment, the pattern 
of recoiling atoms by absorption image method reflects 
the atomic momentum spectroscopy or momentum dis- 
tribution 0, S 0, H H, 0, 111 , where the moving atoms 
and the static BEC form a matter wave grating. At the 
same time, the scatteringoptical spectroscopy shows the 
gain process with time [l|. To enhance the scattering 
light signal, the optical cavity is applied in the similar ex- 
perimental setting which usually called collective atomic 
recoil lasing (CARL) [l^, [l^, [l^ , where the atoms are 
forced to maintain in the density grating by the optical 
lattice in the cavity. Different from the above case, here 
we consider the scattering superradiance from an array 
of BECs. 

Optical lattices (OL), created by pairs of off-resonance 
counter-propagating laser beams, offer new opportuni- 
ties to investigate quantum information processing and 
strongly correlated quantum matter [20| . The periodical 
potential in an optical lattice forms an atomic density 
grating, hence to study the superradiance in this array 
the coherence of atom both inner site and between sites 
need to be considered. Therefore superradiance has the 
potential to become a method to detect the coherence 
of atoms in an optical potential. This is different from 
the interaction of light and BEC in an OL trap without 
the atoms' recoiling 2l|, where the inter sites atomic co- 
herence was considered and the inner site coherence is 
neglected. 

To study the superradiance in an optical lattice, there 
are several problems to be considered with regard to the 
theory about the superradiance from BEC [9|. First the 



frequency of optical lattice, usually several kHz^ is much 
larger than that of the magnetic trap (tens to hundreds 
Hz) where the effects of trap is usually neglected for its 
frequency is much smaller than the recoil frequency. Sec- 
ondly, we need to calculate the gain in the special emis- 
sion angle. In the case of a magnetic trap, the atomic 
cloud experiences the maximum gain when the mode is 
along its long axis. And it can also be understood by that 
the direction is selected for the least width of momentum 
to get the maximum gain Q . The long axis and the least 
width of momentum, these two directions are the same 
in a magnetic trap, but not necessarily true in an OL 
trap. Lastly but important, since the light scattering de- 
pends on the coherence of different sites, the interference 
of scattered light results in amplification at some specific 
frequency and suppression at the others. It could pro- 
vide us a method to obtain the information of the atoms 
in OL. The density grating formed by optical lattice and 
the grating formed by moving and static atoms give two 
criterions for the optical amplification. Hence, similar 
to the coherent-enhanced imaging where Raman super- 
radiance is used to pro be the spatial coherence of BEC 
in a magnetic trap |18| . the scattering spectroscopy re- 
flects the cooperative radiation of atoms inner-site and 
inter-sites. 



II. GAIN FOR THE CONDENSATE IN THE 
TRAP 

We consider the model that the atomic cloud is pre- 
pared in an optical lattice, as shown in Fig.[TJ The optical 
lattice is placed along the z-axis with M-sites centered 
at origin. The lattice constant is oq = A/2, the length of 
condensate is L = Map, and the pumping laser incident 
with wave vector ko propagates along its short axis y. 

In order to investigate the superradiant gain, we adi- 
abatically eliminate the excited state for the far-off res- 
onant pump laser and use the rotating wave approxi- 
mation. Therefore the effective Hamiltonian about the 
coupling between the atomic and electromagnetic fields 
can be written as 



H ~ Ha + Hp 



Hi 



(1) 



*Electronic address: xjzhou@pku.cdu.cn 



where the Ha = / d^r>l'^(r)[p2/2m + V(r)]>l'(r) is the 
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FIG. 1: (Color online) The system sketch. For a cigar shaped 
condensate, the optical lattice is formed along x-axis, with M 
sites. The pumping laser is propagating along y-axis. And in 
the X — z plain, the length in the z-direction is bigger than 
that in the x-direction. 



atomic Hamiltonian, for the classical potential of OL is 
V{y) = Vb cos(2:/ao) with the lattice depth Vq- Hp = 
/ d^k?ia'kb^(k)b(k) is the Hamiltonian of photon, and 
the interaction Hamiltonian is 

H^ = j d^kd3r[%(k)*^(r)bt(k)e'(''''-'') "'t'(r) + H.c], 

(2) 

where ^^(r) is atomic field operator, and b^(k) is the 
annihilation operator for a photon in mode k in the frame 
rotating at the pump frequency wo. Here the photon 
energy = c|k| — Wqi f-nd .g(k) is the coupling coefficient 
for scattering between the pump and vacuum modes. The 
interaction between atoms is neglected because it is too 
small in the time scale. 

The atomic field operator can be decomposed to the 
different side modes as 



*(r, t) 



q 



(3) 



the operator Cq refers to the recoiling atom wave function 
which scatters the q mode light, 'i/'o(r) is the ground- 
state wave-function of the condensate with the chemical 
potential When the OL potential is weak enough, 
we could deal this case as in a magnetic trapQ. And 
when it is extremely strong, we are able to approximate 
the potential as the spatial replication of harmonic trap 
V{v) — (T/o/ao)^^' then decompose the operator to the 
eigenstate of the trap that 



and we define the coefficient 



(4) 



(5) 



ipn is the eigenfunction of n-th level of the trap, and 
the P(n, A) is the Poission distribution with parameter 



A = {q^Jh/2niu)T)'^ = uJrjiiJTi here ujr = hq^ /2m is 
the recoiling frequency. Considering the first-order side 
modes, then the interaction Hamiltonian becomes 

H. = Y.Y.j d'k[?i.9(k)pq(k)A„ct6t(fc)co+-ff.c.], (6) 



q/O n^O 



with 



Pq(k) = j d3r|,/,o(r)|2 exp[-i(k - ko + q) • r] 



(7) 



is the Fourier transform of the ground state density dis- 
tribution centered at kg — q. 

As the experiment shows that there are just several 
side modes dominate the whole scattering process, and 
in order to simplify this problem, we just take one side 
mode q into consideration. Under the Born-Markov ap- 
proximation, the optical field could be obtained that 



h{t) = 6(0)e 



^5(k)Pq(k)A„ctco<5(c.fe). (8) 



Inserting ([8]) into the dynamic equation of atomic field, 
c„ = [H, Cn]/ifi^ we obtain its evolution equation 



dt 



A, 



Gq 4 Co 



N ^ 



(9) 



with the EEC's gain with the q-th mode as 



Gci = N^J d^k\p^{k)\H{\k\ - ko), (10) 

where we assumed that (/(k) is isotropic in the k^ — kz 
plane. The first term on the r.h.s. of Eq.® is the gain 
from the condensate. The second term is the quantum 
fluctuation, which has been discussed in and does not 
affect the superradiant behavior in long time, and we just 
take it as an initial seed. The last term on the r.h.s. of 
Eq. ([9]) is the energy term. 

In order to compare the two cases of pumping the con- 
densate with and without the external potential, we need 
to discuss equation ([9]) by the mean-field approximation, 
replacing the field operator c„ by a c-number c„. By 
transformation c„ — c„ exp(— ia;„t), the equation ^ be- 
comes 



dt 



A, 



Cq Cgcp 

' 2 N 



E 



,^771 6 



i{n — m)ujTi 



(11) 



For the case that potential is turned off when pump- 
ing the condensate by the laser, we assume that there 
is just only one level n satisfying the condition An = 1, 



with eigenenergy uJn = p /2m 



and for the else 



An = 0. Thus on the r.h.s. of the equation pT|) the 
phase factor g»("-™)'^Tt jg always equal to one. And this 
form is consistent with the equation (8) in Q. The gain 
of atomic number is therefore Gq. For the case that 
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the external potential exists, the parameters An are cen- 
tered at uj^/^t and have a standard deviation y^ujj./ujT, 
and we just need to consider these 2y^uJr/LL'T trap levels 
that make A,, ^ 0. And the phase factor e*^"-")'^^* is 
different for different levels. These factors could be ap- 
proximated by calculating the average phase difference 
ujT\/^r/^Tdt = ^LOTijJrdt in a during time dt. Thus the 
different phases in different levels lead to that the sum 
in the r.h.s. of equation PT|) is smaller than the case 
that they are the same phase, in the previous case. This 
results in the gain loss, given by 



Gq — Gq OC — 



(12) 



A typical value of the gain is Gq = 4 x 10^ for 
/ = lOOmW/cm'^ and N = 10^, which is close to the 
frequency of the atom in the optical lattice trap. Thus 
when the ujt is large enough to dephase the coherence of 
condensate and side mode, the gain of light is suppressed. 
Since the effect of the trap is just a shift in the gain, in 
the below we mainly discuss the gain without the trap 
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FIG. 2: (Color online) The sampling and the gaussian factor. 
The solid line is the sampling function for M=10, which is the 
result of interference between different site, and the gaussian 
factor is the result of single site amplification of light, which 
are drawn in dash line for = O.lao, and in dash dotted line 
for = ao- 



III. GAIN FROM AN ARRAY OF 
CONDENSATES IN RELEASED TRAP 

Now we consider the case that the pump beam im- 
mediately incidents after switching off the potential, so 
the external trapping potential could be neglected. The 
wave function of ground state for the i-th site is the Wan- 
nier function Wi{r), which is approximated by the gaus- 
sian function cxp(— ^^^-^ r|/2(Tj) with the half width 
of the wave function aj in j direction {j = x,y,z). 
Thus the ground-state wave-function in the optical lat- 
tice is given by ipoir) = Cnor^Wi{x,y, z), where the 
Cnor = (/ d'^r] X]i=i Wi{Y)\'^Y/'^ is the normalization fac- 
tor. For one site, we assume <Jx, <^y, the maximum 
gain is the z direction because the photon could expe- 
rience the most atomic amplification Here, we need 
to carefully discuss the gain for the whole atomic cloud. 
Since the M sites are placed along x-axis and centered 
at origin with lattice spacing ao, then the atomic density 
can be expressed as 

M M-1 

\Mr)\' = ClAT. + 2 w,{r)w,+^{r)], (13) 

i=l i=l 

where the first term on the r.h.s. describes the atomic 
density of site i. The second term is the overlapping 
between neighboring sites, which is considered only when 
the wavefunction of one site is wide enough to overlap its 
neighbors. The Fourier transformation of the density at 
q = ko is 

(14) 



We proceed to calculate the superradiant gain in pO|) . 
The factor in Pq(k), sin(-^^^fca;)/ sin(^A:2;), as shown 
in solid line of Fig. [2] gives the density profile a sam- 
pling, which means that the Pq(k) is not zero in re- 
gions with half width 211 /Mao and these regions sep- 
arating 27r/ao from each other. Moreover, the factor 

of Pq(k), exp(— ^4^), is centered at ko — q and 

spreads in ki direction, which is significant in the region 
\ki — ko-\- q\ < l/ci. The gaussian factor of Pq(k) is the 
dashed and dot-dashed line in Fig. ^ which shows that it 
is wider when the width in the single site wave function 
is smaller. We assume that ko 3> l/cj, 3> I/oq, which 
means that the atomic momentum is narrow enough that 
all its components could contribute to the optical amplifi- 
cation iOj. And in the region of Pq(k) where the gaussian 
factor is significant, we can approximate the surface of 
the sphere |fc| = |A:o| as a plane tangent to the sphere, 
and the integral in (fTU]) is equivalent to be the integral 
on this plane. So the integral when q = kg + \kQ\9, where 
9 = cos 0kx + sin Oh^, , is a unit direction in the — k^ 
plane. 

Here we mainly consider the maximum gain of light 
in two extreme directions, 9 = corresponding to the z 
direction, and = 7r/2 for a; direction. For 9 = 7r/2, the 
gain is 

G, = Go , (15) 

where 

Go = fjGi,-[l + exp(-^)], (16) 
is related to the normalization factor Cnor, the width 
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in y-direction ay and the factor [1 + exp(— ^r)] which 
reflects the coherence of neighboring sites. Eq. ([T5l) shows 
that the gain Gx is proportional to the number of site 
squared , which is the result of cooperative radiation. 

Because the non-zero region of i5-function in (|10|) is 
a plane parallel to the kx — ky plane when 6* = 0, the 
profile width of Pq{k) will affect the result of Gq. Hence, 
we need to discuss two condition ax < clq and ax > ao- 
For the case that ax < ao, Pq(k) has ao/ax side bands, 
after summing all these side bands whose half width is 
1/Mao, we get the superradiant gain given by 



= GoM 



2 00 1 ^ M 

— Ltq ■ 

ax Mao Gx 



(17) 



For the case that ax > ao, there is just one non-zero 
region of /9q(k), thus, the superradiant gain is 



G, = GoM' 



1 r< 

AI qq ao 



(18) 



In both cases, the gain Gz is proportional to Af , due to 
the incoherent sum of different sites. 

When ax < ao, the gain ratio for the two extreme 
direction is that ^ = -^^^ = L/a^ which is the aspect 
ratio, consistent with the theory without OL trap. When 
ax > do, the gain ratio becomes ^ = ^^^^ which is the 
effective length ratio in this two directions. It should be 
noted that this theory is sound under the condition that 
ko 3> 1/ax- A typical value is kg = 27r/780nm, hence 
ax S> 780nm. If we need the gain in z direction larger 
than that in x direction, we need az ^ M-780nm, which 
is hard to be realized in experiment. Thus the radiation 
usually takes place in the x direction. 



IV. SPECTROSCOPY OF SUPERRADIANT 
SCATTERING 

In the previous section, we understand that the gain of 
light is usually propagating along the a;-axis. Considering 
q = 5X in the x-direction, for the different q the gain can 
be expressed as 



Gq = — exp[ 

a;.. 



a2(fco+g)^sin2^^is4H±£) 



2 ao(fco+g) 



(19) 



sm 



In this equation, we know that the maximum gains 
emerge at = mr. In other words, the gain has 

maximum around fco + 2mr/aQ with separation 27r/ao. 

The spectroscopy of different width of single site wave 
function is plotted in Fig. [3l As shown in subfigure (a), 
when ax > ao, there is only one peak in the gain. Sub- 
figure (b) shows that when the ct^: ^ ao, there are side 
bands. The radiant light has side bands when the wave 
functions of neighboring sites are not overlapped. The 
reason is that atoms in different sites are pumping by 
the same phase light and become to the same phase 
dipolc. The radiant light which is propagating along 




FIG. 3: (Color online) The gain with a factor a^/Go, in x 
direction for the ax ~ ao (a) and ax = O.lao (b). The gain has 
satellites when the width of single site wave function is large. 
When it is narrow, the satellites disappear. Here, M — 10. 




FIG. 4: (Color online) The ratio of maximum gain to second 
maximum gain in x direction versus the one site wave packet 
width. Here, M = 10. 



the lattice has a phase difference in neighboring sites, 
(k — ko) • aox. Thus radiant lights with different fre- 
quencies will have different gains by the averaging over 
the whole lattice. The constructive interference will sin- 
gle out the frequency component satisfying the condition 
(k — ko) • aox = 2mr to amplify, and other components 
will be suppressed due to the destructive interference of 
M sites. 

For the larger width of single site wave function, the 
gain of side bands is smaller. Thus by measuring the 
side band gain could give us a method to obtain the in- 
formation about the width ax- Fig. 3] shows the ratio 
of maximum gain to the second maximum gain versus 
ax- By the spectroscopical measuring, we could obtain 
the information of the width of wave packet, which is 
relevant to the potential quantum phase transition. 
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V. DISCUSSION AND CONCLUSIONS 

111 the BEC superradiant experiment a photon is 
scattered by one atom in the BEC, and this atom acquires 
the recoil momentum. The moving atoms and the static 
BEC form a matter wave grating, which enhances the 
same direction scattering. Due to the mode competition, 
the highly directional emissions of light travel along the 
long axis of the condensate. The stability of the relative 
phase between different atomic matter waves determines 
the coherence time of the matter wave. On the other 
hand, in the coherent atomic recoil lasing (CARL), the 
situation is changed by the presence of the cavity, the 
coherence is preserved as the relative phase of the cavity 
light-fields, and it is independent of the atomic motion 
while given by the cavity linewidth. 

Different to these experimental scheme, here we ex- 
tend the theory of superradiance of BEC [10[ to the case 
of in OL trap. In this trap, an array of atoms from the 
optical lattice form a density grating, and the superradi- 
ance gain are calculated in the quantum theory. In this 
theory, we consider inner-site and inter-site coherence of 
atoms. Only the scattering light satisfy the condition 
that (k — kg) • aqx = 2m: will be singled out to be am- 
plified and the other components will be suppressed in 
different extent. Together with the grating formed by 
the static and moving condensate, both gratings give fre- 
quency selection rules. It is similar to a diode laser with 
internal and external cavities. Only the light is resonant 
to both cavities will be amplified. 

The motion of recoiling atom in the high-frequency OL 
trap will dissipate the coherence resulting in a loss of the 
optical gain proportional to ^Jlot^t- In the magnetic 



trap, the trap frequency is smaller, and the loss can be 
neglected. It can inhibit the collective radiation when the 
trap frequency is high enough. By calculating the ratio 
of optical gain in the two extreme direction, we show 
that the gain is proportional to the length that the light 
travels in the condensate. 

Depending on the lattice depth, the wave function for 
one site overlap differently with its neighboring sites. 
When the OL potential is low enough, the wave func- 
tions of neighboring sites fully overlap, just like a con- 
densate in the magnetic trap. When the OL potential is 
high enough, the wave functions of neighboring sites are 
separated. The different overlap results in the different 
scattering spectroscopy. Thus the spectroscopy provides 
us with a new method to detect the coherence of differ- 
ent sites. This spectroscopy method offers much more 
precision than the absorption image method. Moreover, 
unlike the time-of-flight method, which is usually used to 
detect quantum phase transition [2^, the spectroscopy 
method is a non-destructive method. More understand- 
ing of this mechanism can help to understand the self- 
organization, especially how the long-range order arise 
in the self-synchronization process. Superradiance may 
be helpful to detect the phase transition between the su- 
perfluid (SF) and Mott-insulate(MI). 
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